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Using a model Hamiltonian, we discuss how we could interpret data obtained from transition 
metal if-edge resonant inelastic x-ray scattering (RIXS) experiments. By analyzing the creation of 
valence excitations from the screening of the core hole and calculating corresponding RIXS spectra 
for metals and insulators, we find that the probability for excitations depends not only on the total 
energy but also on the asymmetric screening dynamics between electrons and holes. 
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I. INTRODUCTION 

The inelastic x-ray scattering cross section from a solid 
state system is enhanced by orders of magnitude, when 
the incoming x-ray energy Wj n is tuned to a resonance 
of the system. For resonant inelastic x-ray scattering 
(RIXS), an electron is excited from a deep- lying core 
state into the valence shell. In particular, the RIXS at 
transition metal if -edges, which is the focus of this work, 
has drawn significant attention recentlyJ^' 3 -'^' 6 -'^^^ 
As illustrated in Fig. [TJ transition metal if -edge RIXS 
involves a dipolar transition from the Is shell into the 
wide Ap band [Fig. QJa)], the screening of the Is core 
hole by the 3d valence electrons [Fig. QJb)], and the re- 
combination of the 4p electron and the Is core hole with 
the emission of a photon [Fig.QJc)]. Unlike RIXS at other 
edges (e.g. transition- metal L-edges), the dipole matrix 
elements in transition metal if -edge RIXS at a particular 
incoming photon energy (u>i n ) and momentum (q in ) are 
constant factors, which makes transition metal if -edge 
RIXS a direct probe for the screening dynamics. The mo- 
mentum (q ou t) and energy (w ou t) dependence of the out- 
going x-ray resulting from the radiative decay of the core 
hole allows measurements of the energy u> = uj- m — u> out 
and the momentum q = qj n — q out of elementary excita- 
tions near the Fermi level or across the gap. The RIXS 
has several advantages compared to other spectroscopic 
probes. In contrast to angle-resolved photoemission, 
RIXS is bulk-sensitive due to the large penetration depth 
of x-rays. The use of a particular resonance makes RIXS 
chemically selective. As opposed to x-ray absorption, 
there is no core hole present in the final state. The RIXS 
has provided unique insights into, e.g., the momentum 
dispersion of charge excitations and magnetic excitations 
in high-T c superconductors and related systems ^di^il 
orbital excitations^^ and magnetic-ordering dependent 
transfer of spectral weight in colossal magnetoresistive 
manganites^ 

Numerical calculations of the RIXS spectrum have 
been done on, e.g., finite clustersj 13 i 14 ' 15 In certain lim- 
its, the RIXS cross section can be related to the dynamic 
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FIG. 1: Illustration of the RIXS process studied in our work, 
that is, transition metal oxide if -edge RIXS. This process is 
distinguished from other types of RIXS, such as transition 
metal oxide L-edge RIXS (Figures 1 and 2 in Ref. Iljfb or 
if -edge RIXS for diamond or graphite (Figure 3 in Ref. 



structure factor S'q(w), either by lowest order perturba- 
tion theorjs 3 -^ or by an approximate representation of the 
intermediate-state dynamics. 16 Within these approxima- 
tions, the probability for the intermediate state core-hole 
potential to create excitations, such as electron-hole (e-h) 
pairs, mainly depends on the total energy of the shake- 
up process. In this work, we show that RIXS depends on 
the energies of the constituents of the excitation, leading 
to a strong asymmetry in electron and hole responses to 
the core hole for the parameter range where the screen- 
ing dynamics become important. We demonstrate that 
the RIXS spectral line shape is determined by the nature 
of the intermediate-state core-hole valence-shell-electron 
excitons. 



II. MODEL HAMILTONIAN AND METHODS 

The total Hamiltonian for the RIXS process shown in 
Fig. [1] has electronic, photonic, and electron-photon in- 
teraction parts, 



iitotal — iicl + H p hi + ii i_pht- 
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The RIXS intensity is obtained from 



/(Win, Wout) = -r- ^2 l^'^/l 2 5 (^/ + W ° ut _ ^ _ ' 
/ 

,(2) 

where Ui— ►/ is the transition amplitude from the initial 
state \i) to a final state |/). In terms of the initial 
many-electron energy Ei, many-electron eigenstate |n) 
and eigenenergy E n of H e \, and the FWHM intermediate- 
state lifetime-broadening 2r, the transition amplitude is 
given by 
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The Coulomb interaction between localized Is holes 
and 3c? electrons, which is the central mechanism for tran- 
sition metal fiT-edge RIXS, is U c = 4 - 7 eV (Refs. 0, Q3, 
[H, and HI). The excited Ap electron is highly delocalized, 
and its Coulomb interaction with Is holes or 3c? electrons 
is neglected as an approximation, leading to 



H c \ — Hi s + H 3 d + -ffis_ 



3d ' 



H 



4p. 



(4) 



where Htg, H^d, and H^ p are Hamiltonians for Is, id and 
Ap levels, and Hi s -3d represents ls-3c? Coulomb interac- 
tion. The Ap level Hamiltonian H± p is decoupled from 
the other terms, and can be replaced with a constant £4 P 
for the intermediate state, which we drop from H c \. At 
a particular resonant ui- m for the if -edge RIXS, the ma- 
trix element of iJ e i_ p ht ~ p ■ A between Is and Ap is a 
constant factor, leading to 
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where s} is a Is hole creation operator. Since transition 
metal K-edge RIXS probes initial and final states with- 
out core holes and intermediate states with core holes, \i) 
and |/) are the energy eigenstates of H e \ with s}s = 0, 
that is, i?3d, whereas \n) are the energy eigenstates of H e \ 
with s}s = 1 . Equation {5]) is evaluated by expanding 3d 
parts of \n) in the basis of the energy eigenstates of H 3 d, 
which include 3d parts of \i) and |/). 

In terms of a 3c? electron creation operator d\ the 3d 
Hamiltonian for metals and insulators is 



H. 



a/i 



k 



£kd\dk, 



(6) 



for which we consider a constant DOS without and with 
a gap (A) for metals and insulators, respectively. The 
Hamiltonian includes ./V one-electron 3c? levels. The num- 
ber of 3c? electrons is L. For independent electrons with 
a single orbital degree of freedom at each site, the two 
spin channels are equivalent, and two bound states are 
present at the site with the core hole. However, for 



many transition-metal compounds, the Coulomb inter- 
action between valence electrons is often non-negligible, 
resulting in the Mott insulators. For these Mott insula- 
tors, we consider the Hubbard Hamiltonian, 



rrHubbard 
H 3d 



e h d\ a d Ka + U 4r c Mlj, < % ( 7 ) 



which includes the effect of the on-site 3c?-3c? Coulomb in- 
teraction. Spin and site indices are a and i, respectively. 

Because the Is band is much narrower than the 3g! 
band, the total RIXS intensity is, in a good approxi- 
mation, proportional to the intensity calculated for a Is 
hole created at, for example, site i = (see footnote 
Therefore, we use Is hole Hamiltonian, 



His = SlsS's 



(8) 



The following identity shows that the on-site Coulomb 
interaction with a localized core hole scatters any 3c? state 
with a wavevector k into the same or any other 3c? state 
with a wavevector k' with an equal probability, 



U c d\ =Q d i=0 sl =0 ± 



k,k' 



(9) 



For our model, this corresponds to the scattering from 
any 3c? level into the same or any other 3c? level with 
the same strength of —U c /N, which is expressed in the 
Hamiltonian, 



H 



= -J2l§ d k d k'S 



ls-3d 



k,k' 



N 



(10) 



Our model Hamiltonian for metals and insulators with 
the core hole, that is, H c \ = H 1 ^/ 1 + His~3d + -His. 
with s}s — 1 , is an exactly-solvable independent-particle 
Hamiltonian, and has been solved in the context of x- 
ray edge singularity 20 ' 21 The intermediate /-electron 
eigenstates and eigenenergies are obtained by diago- 
nalizing H c \ exactly and filling L one-electron eigen- 
states. To calculate the RIXS intensity, we classify final 
states according to the number of e-h pairs, for exam- 
ple, the Fermi sea (i.e., the initial state \i)), states with 
one e-h pair d k dk< and states with two e-h pairs 
d k d k' ^fc<^fc' where fc< and fc> represent 3c? electron 
levels below and above the Fermi energy, respectively. 
We examine the contribution of each class of states to the 
RIXS intensity, as we increase the number of e-h pairs in 
the final states. 

For the Mott insulators, the Hubbard model cannot be 
solved exactly. However, it is known that the main fea- 
tures of the Mott insulators are understood by following 
Hubbard's own approximation 22 for i/^ ubbard , which we 
extend to calculate the RIXS spectral line shape for the 
half-filled Mott insulators. First, we project the opera- 
tor for a particular spin onto different occupations of the 
opposite spin, 



4a = (! - m ,-a) d \o + n %-°4o 



(11) 
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The first operator (singlon) g] a = (l—ni-<j)<V ia describes 
the motion of the electron in the absence of the electron 
with the opposite spin. The second operator (doublon) 
b\ a = n^- a d\ a describes the motion of the electron when 
the site is doubly occupied. The local potential in the 
presence of the core hole at i = is written as 

H M -3d = - Uc (4=0 lO -&=0,<7 + &Lo, ct ^=o,<t) • (12) 

<T 

In the Hubbard-I approximation^ the projection oper- 
ators, 1 — Tif,— cr and riij—o-, are replaced by their expecta- 
tion values, resulting in an effective independent-particle 
Hamiltonian. The one-particle eigenvectors for the lower 
and upper Hubbard bands in the absence of the core-hole 
potential are written as 

kkicr) = cos# k |sw) + sin^ k 1 6 k(r ), (13) 
\evva) = sin0 k |sw) - cos6> k |6 k(T }, (14) 

respectively, where k is the wavevector. Although the 
core-hole potential in Eq. fT2")) does not couple singlon 
and doublon states, scattering between the lower and up- 
per Hubbard bands still occurs due to the mixing of the 
singlon and doublon states through electron hopping. In 
the initial state, only the states in the lower Hubbard 
band are occupied, 

K)=ni<w. ( i5 ) 

kcr 

For the intermediate states, the effective independent- 
particle Hamiltonian is solved in the presence of the core- 
hole potential, leading to one-particle eigenstates |e mCT ), 
where the index to = 1, 2, 3, • • • labels the eigenstates in 
the order of increasing energy. We note that the wavevec- 
tor k is no longer a good quantum number and that the 
lower and upper Hubbard bands are no longer well de- 
fined because of the presence of a bound state at the site 
with the core hole. Many-electron intermediate eigen- 
states and energies are obtained by filling an appropriate 
number of the states, \e ma ). By following the same pro- 
cedure as for metals and insulators, we find the RIXS 
spectral shape for the Mott insulators. 

III. RESULTS 

To find RIXS spectrum, all possible intermediate states 
should be considered in principle, as done for small clus- 
ters in Refs. [IH and0- However, for most RIXS exper- 
iments, the incoming x-ray energy is tuned at the ab- 
sorption threshold for the resonance enhancement effect. 
Therefore, without the short lifetime of the intermediate 
state and the energy uncertainty, the intermediate state 
would be purely the lowest energy state, in which the lo- 
calized core hole is screened by the 3d electrons. To find 
the contribution of different intermediate states to the 
RIXS intensity, we calculate a part of the RIXS matrix 
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FIG. 2: X-ray absorption spectroscopy (XAS) spectrum for a 
metallic DOS with a bandwidth W = 2 eV and a core-hole 
potential U c = 5 eV at half-filling. 

element (fT-js^i) in Eq.([5]), namely the projection of the 
intermediate states onto the state of the Fermi sea with 
a core hole, which corresponds to the x-ray absorption 
intensity, 

/xasM oc \(n\s}\i)\ 2 S(u -E n + Ei). (16) 

n 

The x-ray absorption intensity, calculated for a square 
density of states with a band width W = 2 eV and U c = 5 
eV at half-filling, is displayed in Fig. [21 which shows two 
distinct spectral bands corresponding to the XAS final 
states with an almost full (screened) and an almost empty 
(unscreened) bound state, separated by approximately 
U c (Ref. I20I ). Experimentally, the energy difference be- 
tween the onsets of resonant inelastic x-ray scattering for 
a screened and an unscreened core hole is 5-7 eV (Ref.!), 
which corresponds to U c . For the lower spectral band, the 
XAS spectral weight is dominated by the transitions into 
the lowest energy intermediate state, |ni ow ), whereas the 
higher-lying states in this band have little component 
of the Fermi sea with a core hole, s}\i). It shows that 
the lowest energy screening state, indicated by an arrow 
in Fig. [21 dominates the matrix element. The smaller 
satellite peak at higher energy corresponds to the excita- 
tion of the unscreened state. Since the core-hole lifetime 
broadening T is significantly smaller than the energy dif- 
ference between the screened and unscreened absorption 
states, the probability for the unscreened state excitation 
would be low. Even if the unscreened state is partially 
excited, its matrix element for RIXS is smaller than that 
for the screened state. Based on these arguments, we 
propose that the screened state is the dominant inter- 
mediate state contributing to the RIXS intensity, and we 
replace the sum over the intermediate states with a single 
lowest energy intermediate state, ni ow : 

I(Ato) oc Y K/kKw}| 2 S[Au; - (Ef - Ei)]. (17) 

For metals and insulators, using indices k and to, which 
are reserved for 3d-level one-electron eigenstates with- 
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out and with the core hole respectively throughout this 
paper, we express operators for the intermediate one- 
particle eigenstates with energy e m (m = 1,2, ...,N in 
the order of increasing energy) as 

Di = Y. a ^ d i ( l8 ) 

k 

The numerically exact lowest energy intermediate state 
l^iow) is found by filling up the lowest L one-electron 
energy eigenstates of H c \ with s}s = 1 . 

The RIXS intensity in Eq. [T7] calculated for the final 
states with up to three e-h pairs for a half-filled metal- 
lic DOS is shown by the upper (black) solid line marked 
with A = in Fig. Ufa). The full bandwidth and the 
core-hole potential are chosen as W = 2 eV and U c = 5 
eV, respectively. The RIXS intensity calculated only for 
the final states with two e-h pairs is shown by the lower 
(green) solid line marked with A = in Fig. |3fa), which 
is small compared to the total intensity. The contribution 
of the final states with three e-h pairs to the RIXS inten- 
sity is negligible. The results show that the final states 
with one e-h pair contribute dominantly to the RIXS in- 
tensity, which means that the lowest intermediate state 
is approximately a combination of the Fermi sea and the 
states with one e-h pair, 

t K>+ E Aow^AkML'^H)- (is) 

k > ,k < 

The first term gives rise to a large (ni ow |s^|i}, the XAS 
part in the RIXS matrix element in Eq. (O as discussed 
earlier. The remaining terms with one e-h pair give a 
large (/|s|ni ow ), the other part in the RIXS matrix ele- 
ment. 

The electron and the hole states in \n\ ow ) are entangled 
in general. We examine whether unentangled electron 
and hole states serve as an approximation, that is, 

|niow) ~ atum<L\i) + 

5>fow(fc>)4> EAow(fc<H< a t l<>.(20) 
fc> / \fe< / 

We calculate the electron and hole distribution functions 
for the lowest energy intermediate state using a m k de- 
fined in Eq. HID, 

n t>( £ k>) = E i a ™ fc >i 2 ' ( 21 ) 

m— 1,L 

= E l fl ^<! 2 ' ( 22 ) 

m=L+l,N 

which are shown in Fig. Efb) for A = 0. If Eq. ^ is 
exact, we would have 

n e k> (s k> ) cx |A e ow(fc>)| 2 , (23) 
nt(e k< ) cx K w (k<)\ 2 (24) 




FIG. 3: (color online) (a) RIXS spectra calculated for a square 
density of states with a band width W = 2 eV and U c = 5 
eV at half-filling for a metal [upper (black) solid line for A = 
0] and insulators (A > 0). For the insulators, a gap (A) 
is opened at the center of the band, and the occupied and 
unoccupied bands have the same width W/2. For the metal 
(A = 0), the lower (green) solid line shows the contribution of 
the double electron-hole pairs to the RIXS spectrum and the 
upper (red) dotted line gives the approximate spectrum using 
Eq. (|25|l . The lower (green) dotted line is a shifted energy 
axis, (b) Distribution functions of holes (n k< ) and electrons 
(n| ) below and above the Fermi level for different values of 
A, created by the intermediate-state core-hole potential. 

and 

I(Au}) = A E nU^X^H^-^+^I 25 ) 

k> ,k< 

Therefore, the comparison of the RIXS intensity calcu- 
lated from Eq. (|2"5|) with the one calculated from Eq. (jTTJ) 
reveals whether Eq. (|20| is a good approximation. The 
RIXS spectrum calculated from Eq. |25|) is shown with 
the upper (red) dotted line for A = in Fig.[3ja), which 
shows a good agreement with the upper (black) solid line, 
indicating that a separable e-h pair state is a good ap- 
proximation for the lowest energy intermediate state. 

We examine the characteristics of the electron and hole 
distribution functions, n e (e/ £> ) and n' l (£fe < ), shown in 
Fig. OJb) for a metal (A = 0). First, the distribution 
functions show a strong singularity near the Fermi level, 
which is responsible for a similar singularity in the RIXS 
intensity, I(Aoj). These features correspond to the low 
energy e-h excitations right near the Fermi energy, re- 
lated to the x-ray edge singularity^ However, the elas- 
tic peak in RIXS overlaps with this edge singularity, and 
higher energy features are more important. Away from 
the Fermi energy, hole spectral weight, n k (£fe<), dimin- 
ishes very rapidly, whereas the electron spectral weight, 
n| (£&>), remains nearly constant. We interpret this 
difference in the following way. In the lowest energy in- 
termediate state, the electron screens the localized core 
hole, meaning that the excited electron should be local- 
ized in the real space. That corresponds to the state with 
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a nearly constant amplitude in fc-space, or a constant am- 
plitude in the energy space since we assume a constant 
electron DOS. The excited 3d hole does not screen the 
core hole, but is in a delocalized state right below the 
Fermi energy to minimize its kinetic energy, which cor- 
responds to the spectral weight concentrated near the 
Fermi energy. Such asymmetry gives rise to the higher 
energy features in RIXS. The RIXS intensity for A = 
in Fig. [3] is large up to 1 eV, i.e., the width of the un- 
occupied band, beyond which the intensity quickly de- 
creases to zero. This asymmetry differs significantly from 
the existing approaches in Refs. [1, @ andll6l which pre- 
dict that the RIXS spectral line shape is proportional 
to the dynamic structure factor S(u) (i.e., the convolu- 
tion of the occupied and the unoccupied DOS) weighted 
by the energy of the shake-up structure. According to 
these studies, the RIXS for a square DOS would be given 
by S(u>), which has a triangular shape, multiplied by a 
U 2 /uj 2 (Refs.iandi) or a U 2 /(U c -ui) 2 (Ref.M) energy 
dependence. However, Fig. (3Ja) clearly shows strong de- 
viations from a simple triangular line shape. The results 
in Refs. and0 are correct in the limits of U c — -> and 
U c — > oo, respectively. For intermediate U c , the screen- 
ing dynamics becomes more important, leading to the 
asymmetry in electron and hole excitations. 

We also carry out similar calculations for insulators. 
The results are shown in Fig. [3] for A > 0. For a small 
gap, the singularity disappears and the RIXS spectral 
line shape resembles more closely the unoccupied DOS 
than S(lu). For larger gap values of the order of the 
bandwidth, the spectral line shape develops into a tri- 
angular shape. The electron-hole asymmetry between 
nj! and becomes weaker, but persists. It shows 
that the widths of RIXS peaks for metals and insulators 
with gaps smaller than the band width roughly corre- 
spond to the widths of unoccupied bands. Therefore, 
RIXS probes mainly unoccupied bands, complementing 
ARPES, which probes occupied bands. 

The RIXS for the half-filled Mott insulators is calcu- 
lated by the same procedure used for metals and insu- 
lators. For the lowest intermediate state, the lowest N 
one-particle eigenstates are filled, 



|niow) = Yl \ £ mo) 

cr=t,J.,m=l,— ,N/2 



(26) 



De-excitations from the lowest intermediate state are 
obtained by representing |e mcr ) in terms of |£k£cx) and 
|£k{7cr)j and calculating the overlap with final e-h pair 
excitations, such as 



I/) = Is^Ua') J| kkicr). 

ko"7^k< a' 



(27) 



The RIXS spectrum is calculated according to Eq. (fTT)) . 
The results for U = 1, 2, and 3 eV are shown in Fig.[4ja), 
where we consider final states with up to two e-h pairs. 
Electron and hole distribution functions for the Mott in- 
sulators in Fig. HJb) show even stronger asymmetry than 
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FIG. 4: (color online) (a) RIXS spectra calculated for the half- 
filled Mott insulators with the on-site Coulomb interaction 
U = 1,2, and 3 eV (U c = 5 eV, W = 2 eV). (b) Distribution 
functions of holes and electrons created by the intermediate 
core-hole potential for the Mott insulator with (7 = 1 eV. 



those for metals and insulators. The width of the RIXS 
peak is about W/2, which is the width of the unoccupied 
(upper Hubbard) band. More discussions are given in 
the next section. 



IV. DISCUSSIONS 

We analyze the relations between the one-particle 
eigenstates with and without the core hole to understand 
the results obtained in Sect. III. As mentioned in Sect. II, 
the Hamiltonians with the core-hole potential have been 
studied in the context of x-ray absorption edge singular- 
ity ) 20 i 21 and it is known that the coefficient relating the 



two basis sets in Eq. (1181) i s 



£fc — En 



(28) 



where C m is a normalization factor and e m 's are the so- 
lutions of 



1 1 



^ Ft. — 



U c N ^-- J e k — e, 

K 



(29) 



Figure 1 in Ref. [2CJ illustrates how e m is determined 
graphically for a given set of e^s and a value of U c . 

First, we analyze the results for metals. For a metallic 
DOS with N states separated by an equal, small energy 
difference, the effect of the strong attractive Coulomb 
potential — U c is to pull out a bound state below the 
bottom of the band and shift the energies of the remain- 
ing N - 1 states (see Fig. 1 in Ref. HO). The calcu- 
lated DOS in the presence of the core-hole potential is 
shown with the gray (green) solid line in Fig. \5j[b). The 

bound state energy is s m —\ = —U c — j^j-- The en- 
ergies for the rest of the states are e m = — 
where the indices m and k have the same value and the 
phase shift is S m = n, . . . , for Ek = — . . . , ^- in the 



large U c and iV limit. As indicated by Eqs. (|2Tj) and 
(|22p . a m fe > 's with m = 1, L correspond to electron ex- 
citations, because they represent levels unoccupied ini- 
tially but occupied in the intermediate state. Similarly, 
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amfc<'s with m = L + l,...,iV correspond to hole exci- 
tations. The expression of a m ^ — C ot /(e/- — e m ) indi- 
cates that the squared amplitude of a state A; in a state 
to decreases as their energy difference increases, propor- 
tional to (ek — £ m )~ 2 - Therefore, if the intermediate 
one-electron energy e m is well separated from a band 
of £fc, like the bound state to = 1, all states in the 
band would contribute almost equally to that interme- 
diate state. In contrast, if the intermediate state en- 
ergy e m is right next to some e^'s, like the intermediate 
states in the continuum band, those nearby Sk states con- 
tribute dominantly to that intermediate state. Therefore, 
as indicated by the long (blue) arrow between Figs. [3a) 
and [5jb) and the dotted line in Fig. EIJa), the interme- 
diate bound state includes electron excitations with an 
almost constant distribution function above the Fermi 
energy. The intermediate state in the continuum in- 
cludes electron or hole excitation only if the intermediate 
state is right near the Fermi energy, as indicated by the 
two short (pale blue and red) arrows between Figs. [5](a) 
and[5lb), giving rise to singular features. This idea is fur- 
ther demonstrated by the (red) dashed line in Fig.[5Jb), 



which represents n e m — J2k> a mk> f° r m ~ 1> ■••jA the 
levels with energies below zero, and = J2k< a mk< f° r 
to = L + 1, JV, the levels with energies above zero. 
The distribution functions n e m and represent how dif- 
ferent intermediate eigenstates contribute to the the e-h 
excitations, whereas ra£ and n\ represent how differ- 
ent initial and final eigenstates contribute to the the e-h 
excitations. 

We apply a similar analysis to insulators, as shown 
in Figs. 03c) and[5](d). The intermediate state DOS is 
shown in gray (green) solid line in Fig. [5^d). It indicates 
two bound states, one below the bottom of the lower 
band and the other within the gap. Other L — 1 and 
TV — L — 1 states are within the lower and upper contin- 
uum bands, respectively, with energies shifted in a similar 
way as metals. Therefore, the L electrons fill right to the 
top of the lower band. In the limit of U c , A 3> W, the 
two bound state energies are — % ± lyt^c + ^ur- 
ther approximation in the limit of U c 3> A gives rise to 
the bound state energy below the bottom of the band, 
e m= i = —U c — A 2 /(4[/ c ), and the bound state energy in- 
side the gap, £ m =L+i — A 2 /(4(7 C ). The (red) dashed 
line in Fig. [Hd) shows n e m for to = 1, ,..,L and 
for to = L + 1,...,N. Because of the gap, the inter- 
mediate continuum states do not include e-h excitation, 
and therefore, the singularity in n% or n\ < disappears. 
The states with substantial electron and hole excitations 
are the bound states, e m =i and e m= L+i, respectively, as 
shown with the arrows between Figs.[5](c) and[5£d). The 
shape of electron and hole distribution functions, n| 
and n\ , are determined by the distance between the in- 
termediate bound states and the initially occupied or un- 
occupied bands. The bound state below the lower band 
is well separated from the initially unoccupied band, giv- 
ing rise to an electron distribution function close to a 
constant. In contrast, the bound state within the gap 
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FIG. 5: (Color online) (a) and (c) show the distribution func- 
tions of holes and electrons in the intermediate state for a 
metal and an insulator (A = 0.3 eV) respectively, similar to 
the ones shown in Fig. [3jb). (b) and (d) show the intermedi- 
ate DOS with gray (green) solid lines and and nj^, denned 
in the text, with (red) dashed lines. Arrows indicate how the 
two plots are related, as explained in the text. For illustrative 
purpose, a smaller core-hole potential (U c = 2 eV) is chosen. 
Other parameter values are identical to those for Fig. 



is close to the top of the initially occupied band, which 
gives rise to a hole distribution function peaked near the 
gap. This difference results in the strong asymmetry in 
the electron and hole distributions. 

The results for the half-filled Mott insulators shown 
in Fig. 3] are interpreted in a similar way. The on-site 
Coulomb interaction splits the DOS into the lower and 
upper Hubbard bands (LHB and UHB) with a gap pro- 
portional to U. With two spin degrees of freedom on the 
site with the core hole, two bound states are pulled be- 
low the bottom of the LHB, as shown with the (green) 
solid line in Fig.[Bfb). The lowest bound state at approx- 
imately U c below the LHB consists mainly of the states 
from the LHB. The second bound state, which is about 
U higher in energy, is made of mainly the states from 
the UHB. For the lowest intermediate state, both bound 
states are filled, giving an energy of — 2U C + U, which 
is equivalent to the energy for two electrons at the site 
with a core hole. After the de-excitation, these bound 
electrons scatter to different states. However, the LHB 
is already occupied in the ground state and the lowest 
bound state does not contribute to the inelastic scatter- 
ing, as indicated by n e m and n 1 ^ shown in the (red) dashed 
line in Fig. [H)(b). In contrast, very few doubly occupied 
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FIG. 6: (color online) (a) Occupation of holes and electrons 
for the Mott insulator with U — 1 eV, similar to the one 
shown in Fig. |4^b). For illustrative purpose, a smaller core- 
hole potential (U c = 2 eV) is chosen, (b) Intermediate DOS 
shown in gray (green) solid line and n„ and nj^, defined in 
the text, in (red) dashed line. Arrows indicate how the two 
plots are related, as explained in the text. 



states are present in the ground state and the second 
bound state leads to the scattering to all the states in the 
UHB. After filling the two bound states, N — 2 electrons 
are left to fill N — 1 states in the LHB in the intermediate 
state. Therefore, the lowest intermediate state has a hole 
at the top of the LHB, and the hole scatters only to a nar- 
row region of the occupied DOS, as shown in Figs. EJa) 
and[Sfb). The RLXS spectral line shape therefore mainly 
reflects the unoccupied DOS, as is clear from the RIXS 
spectra in Fig. 0|a). 

The relationship between the RIXS spectrum and the 
dynamic structure factor is one of the interests in recent 
studies. Our results on the forms of and n e k> for 
metals and insulators allow us to find explicit expressions 
of in certain limits, which can be compared with the 
dynamic structure factor S(lu) directly. For example, in 
the limit of U c » A, W, we approximate ni. ~ U~ 



for 4 < 



< A+W 
— 2 ' 



we approximate 

e k< < -f • From these and Eq. (55} , we find the RIXS 
spectrum and its relation to the dynamic structure factor. 
For A < lu < A + ^, we obtain 



IM ~ 5(c)/ 



where S(ui 

~ S(lu)/ 



lu - A. For A 
W + A 



A 



\ < lu < A + W, 



W 



(30) 



(31) 



where S(oj) ~ A + W — lu. This shows that the rela- 
tion between RIXS intensity and the dynamic structure 



factor is not as straight-forward as expected from pre- 
vious studies. As for the integrated spectral weight, we 
find that it approaches to AL(N — L)/N 2 in the limit of 
U c 3> A ^> W, recovering the results found in Ref. HH 

Materials with narrow unoccupied and wide occupied 
bands provide testing grounds for our results, since our 
theory predicts a narrow RIXS spectrum with the hole 
excitations close to the Fermi level, whereas a simple con- 
volution of the unoccupied and occupied states gives a 
wide RIXS spectrum. Recent calculations-- suggest that 
the undoped manganite LaMn03 has such asymmetric 
band widths (about 0.4 eV and 2 eV for the unoccupied 
and the occupied bands, respectively) due to the layered 
A-type orbital ordering. The RIXS spectrum measured 
for this material^! indeed shows a narrow peak with a 
small dispersion (less than 0.5 eV), consistent with our 
theory. It would be also worthwhile to examine the RIXS 
spectrum for materials with prominent features in the 
occupied bands far below the Fermi level, which would 
distinguish our theory from a simple convolution. Wc 
further note that, when Lu- m is tuned to the satellite (un- 
screened) peak in the XAS spectrum, the situation is re- 
versed: a hole is bound to the core hole and an electron 
is excited near the Fermi energy. Such RIXS tuned at 
the satellite peak would probe occupied bands, instead 
of empty bands. 



V. SUMMARY 

In summary, we have shown that the spectral weight 
of charge excitations in transition metal -ftT-edge RIXS, 
where the core-hole potential is larger than the effec- 
tive bandwidth, is not determined by the total energy of 
the shake-up, but results mainly from the nature of the 
intermediate-state core-hole valence-shell-electron exci- 
tons. We have identified the strong difference between 
excited electrons and holes, which can be used to study 
occupied and unoccupied bands separately. The method- 
ology described in this work can be extended to more 
complicated systems, for example, by incorporating re- 
sults from ab initio calculations, allowing a better inter- 
pretation of transition metal if-edge RIXS. 
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